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ABSTRACT 

We present a spectroscopic study of the stellar populations of Low Luminosity 
Active Galactic Nuclei (LLAGN). Our main goal is to determine whether the stars 
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who live in the innermost (100 pc-scale) regions of these galaxies are in some way 
related to the emission line properties, which would imply a link between the 
stellar population and the ionization mechanism. High signal to noise, ground 
based long-slit spectra in the 3500-5500 A interval were collected for 60 galaxies: 
51 LlNERs and LINER/HII Transition Objects, 2 Starburst and 7 non-active 
galaxies. In this paper, first of a series, we (1) describe the sample; (2) present 
the nuclear spectra; (3) characterize the stellar populations of LLAGN by means 
of an empirical comparison with normal galaxies; (4) measure a set of spectral 
indices, including several absorption line equivalent widths and colors indicative 
of stellar populations; (5) correlate the stellar indices with emission hue ratios 
which may distinguish between possible excitation sources for the gas. 

Our main findings are: (1) Few LLAGN have a detectable young (< 10^ 
yr) starburst component, indicating that very massive stars do not contribute 
significantly to the optical continuum. In particular, no features due to Wolf- 
Rayet stars were convincingly detected. (2) High Order Balmer absorption lines 
of HI (HOBLs), on the other hand, are detected in ~ 40% of LLAGN. These 
features, which are strongest in 10^-10^ yr intermediate age stellar populations, 
are accompanied by diluted metal absorption lines and bluer colors than other 
objects in the sample. (3) These intermediate age populations are very common 
(~ 50%) in LLAGN with relatively weak [01] emission ([OI]/Ha < 0.25), but 
rare (~ 10%) in LLAGN with stronger [01]. This is intriguing since LLAGN with 
weak [01] have been previously hypothesized to be "transition objects" in which 
both an AGN and young stars contribute to the emission-line excitation. Massive 
stars, if present, are completely outshone by intermediate age and old stars in 
the optical. This happens in at least a couple of objects where independent 
UV spectroscopy detects young starbursts not seen in the optical. (4) Objects 
with predominantly old stars span the whole range of [01] /Ha values, but (5) 
sources with significant young and/or intermediate age populations are nearly all 
(~ 90%) weak [01] emitters. 

These new findings suggest a link between the stellar populations and the gas 
ionization mechanism. The strong-[OI] objects are most likely true LLAGN, with 
stellar processes being insignificant. However, the weak-[OI] objects may com- 
prise two populations, one where the ionization is dominated by stellar processes 
and another where it is either governed by an AGN or by a more even mixture of 
stellar and AGN processes. Possible stellar sources for the ionization include weak 
starbursts, supernova-remnants and evolved post-starburst populations. These 
scenarios are examined and constrained by means of complementary observations 
and detailed modeling of the stellar populations in forthcoming communications. 
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Subject headings: galaxies: active - galaxies:nuclei - galaxies: stellar populations 
- galaxies:star formation 



1. Introduction 

Low Ionization Nuclear Emission-line Regions, better known as LINERs, were first 
described as a separate class of galaxies by Heckman (1980). These objects have character- 
istically strong low-ionization forbidden lines of [01], [Nil], [SII] and [Oil] that distinguish 
them from both Seyfert and Starburst nuclei. Because of their small luminosities compared 
to Seyferts and quasars, these objects are often called Low Luminosity Active Galactic Nuclei 
(LLAGN). LLAGN are the most common form of nuclear activity in the local universe. In 
the magnitude hmited survey by Ho, Fihppenko & Sargent (1995, 1997, hereafter HFS97), 
they comprise ~ 1/3 of all galaxies. 

What powers LLAGN and how do they fit in the global picture of AGN? Are they all 
truly "dwarf Seyferts" powered by accretion onto nearly-dormant supermassive black holes 
( "dead QSOs" ) or can some of them be explained at least partly in terms of stellar process? 

These questions have been at the forefront of AGN research for over two decades. Nowa- 
days there is no doubt that a substantial fraction of these nuclei contain bona fide AGN. 
Among the most solid evidence in support of this "AGN model" is the fact that ~ 1/5 of 
LLAGN have broad (few thousand km/s) Ha emission (HFS97), signahng the presence of 
a Broad Line Region, a hallmark of AGN. In fact, some of the most compelling evidence 
for the existence of nuclear accretion disks come from LLAGN with double-peaked broad 
Balmer lines (eg, NGC 1097, Storchi-Bergmann etal. 1997; M81, Bower etal. 1996). X-ray 
observations (Terashima etal. 2000; Ho etal. 2001) and the detection of compact nuclear 
radio emission (Falcke et al. 2000; Nagar et al. 2000) also support the AGN model. Sim- 
ilarities between LLAGN and Seyferts were also identified in the properties of their host 
galaxies (Ho, Fihppenko & Sargent 2003). In the theoretical front, it has long been known 
that photoionization by a typical AGN can reproduce the characteristic LINER emission 
fine spectrum (Ferland & Netzer 1983; Halpern & Steiner 1983). These and other evidence 
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(see Barth 2002 and references therein) have contributed to the general understanding that 
many, maybe most LLAGN are bona fide members of the AGN family, the low luminosity 
cousins of Seyferts. 

Despite this general consensus, for nearly as long as LINERs have existed in the hter- 
ature, it has been suspected that they constitute a rather mixed bag of phenomena (e.g., 
Filippenko 1996; Heckman 1996). LINER-like spectra are known to occur in as diverse envi- 
ronments as cluster cooling flows (e.g. Heckman ctal. 1989; Donahue & Voit 1991), starburst 
driven super-winds (Heckman, Armus, & Miley 1990), and the bulges of early-type galaxies 
(e.g. Heckman 1996) besides the nuclei of galaxies. It is therefore likely that a variety of 
ionization mechanisms arc capable of producing a LlNER-type spectrum, and indeed several 
models not involving photoionization by an AGN have been proposed to explain the emis- 
sion hue characteristics of these objects. Among them are shocks (Heckman 1980; Dopita 
& Sutherland 1995), which may be associated with SN remnants, photoionization by un- 
usually hot young stars (Fihppenko & Terlevich 1992; Shields 1992; Barth & Shields 2000), 
or by post-AGB stars (Binette etal. 1994; TanigTichi, Shioya & Murayama 2000). Even in 
nuclei where an AGN is conclusively known to exist, it is unclear what is its role in the 
gas ionization and in the overall energetics, since the AGN may coexist with one or more 
of these other ionizing agents. For instance, CIV and SilV PCygni lines from O stars were 
detected in the space- UV for a few LLAGN by Maoz etal. (1998), lending support to the 
idea that at least some of these objects may be powered by young stars. Yet, in some of 
these same nuclei, like NGC 4569 and NGC 4594, an AGN has been detected by Chandra 
observations (Ho etal. 2001). These nuclei therefore contain both a compact starburst and 
an AGN. NGC 4303, recently studied by Colina ct al. (2002) and Jimenez-Bailon ct al. (2003) 
is another example. This coexistence of a starburst and an AGN is reminiscent of the star- 
burst + Seyfert 2 composites discussed by Heckman etal. (1997), Gonzalez Delgado etal. 
(1998), Storchi-Bcrgmann, Cid Fcrnandcs & Schmitt (1998), Gonzalez Delgado, Heckman & 
Leithcrcr (2001), Cid Fcrnandcs etal. (2001a) and Joguet etal. (2001). 

The question then arises of what are the relative roles of accretion and other processes 
in the energetics of LLAGN, and how to segregate nuclei where accretion is the dominant 
energy source from those whose emission lines are powered mainly by stars or their remnants. 
HFS97 addressed this issue by dividing LLAGN into two subclasses: "pure" or "normal" 
LINERs and LINER/HII "Transition Objects" (hereafter TOs), whose emission hue ratios 
are intermediate between those of LINERs and HII regions. The fundamental distinction 
between LINERs and TOs is the relative strength of the [OI]A6300 hue (e.g. [01] /Ha). 
This collisionally-excited line arises in warm neutral gas, and as such is weak in gas that is 
photoionized by the soft radiation from O stars, since the transition between the external 
cold neutral gas and internal ionized layers is very sharp. In contrast, an AGN produces 
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highly energetic photons that penetrate deeply into (and heat) the neutral gas. Although 
useful, the division of LLAGN onto LINERs and TOs is largely arbitrary, given that there 
is no clear-cut gap between LINERs and HII nuclei in the [01] /Ha ratio. 

The idea when the TO class was introduced was, in the words of HFS97, that "TOs can 
be most naturally explained as normal LINERs whose integrated spectra are diluted or con- 
taminated by neighboring HII regions" (see also Ho, Filippenko & Sargent 1993; Gongalvez, 
Veron & Veron-Cetty 1999). Accordingly, the most popular models for these sources involve 
photoionization by hot stars in a weak nuclear starburst (Shields 1992; Filippenko & Tcr- 
levich 1992; Barth & Shields 2000), although mechanical heating by supcrnovac following a 
short-duration burst (Engelbracht ctal. 1998; Alonso Herrero etal. 2000) and photoioniza- 
tion by evolved "post-starburst" populations (Taniguchi et al. 2000) have also been invoked 
to explain their emission line properties. 

Given that so many different mechanisms are capable of producing a similar emission 
line spectrum, it is evident that emission line data per se will not be able to distinguish 
between different models. 

In this paper and others in this series we examine the stellar populations in a sample of 
LLAGN in search for clues of their nature. A simple motivation for this is that all alternative 
non-AGN models discussed above share a common trait: They all involve stars in one way or 
another, be they young and massive (a starburst), older and not so massive (post-starburst 
models) or massive but dead (SNR shocks). Characterizing the stellar populations of LLAGN 
may thus provide us with the extra information needed to assess the relevance of the proposed 
scenarios. 

With this central goal in mind, we have carried out a spectroscopic survey of LINERs 
and TOs in the 3500-5500 A interval. This spectral range contains several stellar population 
features of interest, including the WR bump at 4680 A (an unambiguous signpost of the 
presence of young, < 10^ yr massive stars), high order Balmer absorption lines (HOBLs) 
of HI (which signal the presence of evolved starbursts, from 10^ up to 10^ yr), the 4000 A 
break, the Balmer jump and several metal lines typical of old populations, such as Call K, 
the G-band and others. A detailed scrutiny of these features by means of stellar population 
analysis provides a direct test of models proposed for LINERs and TOs. This technique can 
be contrasted to the mostly used indirect approach of comparing emission line ratios with 
predictions from photoionization or shock models. 

Our survey is based on the sample of HFS97, but differs from it in two important aspects: 
Wavelength coverage and spatial resolution. Unlike HFS97 we do not reach the red part of 
the spectrum, which includes emission lines which are the basis of spectral classification. On 
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thc other hand, the HFS97 spectra start at 4230 A, whereas our data goes down to 3500 
A. This near-UV edge contains several key stellar population diagnostics (like the HOBLs) 
which have not been systematically studied in previous LLAGN surveys. As shown here, 
results obtained in this spectral range give us a new perspective on the nature of these 
galaxies. 

In this paper we concentrate in the presentation of the data. Companion papers will 
discuss stellar population synthesis analysis (Gonzalez Delgado et al. 2003, hereafter Paper 
II), spatial gradients of spectral properties, morphology, emission line properties and models. 
This paper is organized as follows: In §2 and §3 wc describe the sample, observations and 
reduction process. In §4 the nuclear spectra are presented. A suite of properties has been 
measured from these spectra, such as equivalent widths and colors. These are presented in §5, 
which also compares our measurements with previously published results. This section also 
presents an objective definition of the pseudo continuum used to measure stellar population 
indices in Bica's system. In §6 we examine the most important finding of this investigation: 
The strong connection between stellar population and emission line properties in LLAGN. 
Finally, §7 summarizes our results. 



2. Sample 

We selected galaxies from the Ho et al. (1995, 1997) sample for this study because it is the 
most complete and homogeneous survey of LLAGN available and provides a representative 
population of galaxies in the local universe. HFS97 divided LLAGN into two subclasses: 
"pure" or "normal" LINERs and LINER/HII "Transition Objects" (TOs), whose emission 
line ratios are intermediate between those of LINERs and HII regions. Some Seyfert galaxies 
also match the usual definition of LLAGN, Lna < 10'^° ergs~\ but in this paper "LLAGN" 
is used as a synonimous for LINERs and TOs. HFS97 distinguish LINERs from TOs by 
the value of the [OI]A6300/Hq; emission hne ratio: They define LINERs as nuclei with 
[OI]/Hq; > 0.17, while for TOs 0.08 < [OI]/Hq; < 0.17. Because of its sensitivity to the high 
energy part of the ionizing spectrum, the [01] /Ha ratio is a good indicator of the presence 
of an AGN. 

LINERs and TOs are further distinguished from Seyfert and HII nuclei on the basis of 
other line ratios: [OIII]A5007/H/3 < 3, [NII]A6583/Hq; > 0.6 and [SII]AA6716, 6731/Hq; > 
0.4. According to these criteria, the HFS97 survey contains 94 LINERs and 65 TOs. We 
have added three other galaxies to these lists: NGC 772 (classified as "H/T2:" by HFS97, 
ie, a hybrid between an HII nucleus and a TO) and NGC 2685 (a "S2/T2:") were grouped 
among TOs, while NGC 6951 (classified as a Seyfert 2 by HFS97 and as a Seyfert 2/LINER 
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by Perez etal. 2000) was grouped among LINERs. With these additions the HPS survey 
contains 162 LLAGN. This sample will hereafter be called the "HPS sample". 

We have observed a subset of 28 LINERs and 23 TOs, which corresponds to ~ 1/3 of the 
HPS sample. Por comparison purposes, we also observed 7 normal galaxies and 2 Starburst 
nuclei, which gives a total of 60 galaxies. Table 1 hsts the galaxies in our sample, along with 
several properties extracted from HPS97. As in HFS97, LINERs arc listed as cither L2 (if 
they contain only narrow lines) or LI. 9 (if a broad Ha was detected), and similarly for TOs. 
The L2/T2 and T2/L2 mixed-classes are used when a galaxy looks more like a LINER than 
like a TO and vice-versa, respectively. These ambiguous cases arise due to inconsistencies 
between the classification implied by different line ratios and measurements uncertainties. 
About 13% of the LINERs and TOs in the HPS sample are in these ambiguous categories (and 
many more have uncertain classifications), which shows that the frontier between LINERs 
and TOs is a rather fuzzy one. 

2.1. LINERs, TOs, weeik and strong-[OI] nuclei: cleissification issues 

In this paper we will adopt a slightly different classification criterion for LLAGN. While 
HPS97 divide LINERs from TOs at [01] /Ha = 0.17, we prefer to place this dividing line at 
0.25. To avoid confusion, we shall refer to sources with [OI]/Hq; < 0.25 as "weak-[OI] nuclei" , 
while [OI]/Hq; > 0.25 sources are classed "strong- [01]". Our motivation to introduce these 
definitions is entirely empirical. We shall sec that the [01] /Ha = 0.25 limit represents much 
better the combined distributions of [01] /Ha and stellar population properties in LLAGN. 
Clearly, this choice is largely irrelevant insofar as emission lines are concerned, since weak- 
[01] nuclei still have [01] /Ha values intermediate between strong- [01] and HII nuclei. In 
other words, weak- [01] nuclei are essentially TOs and strong- [01] are all LINERs. 

Of the 162 LLAGN in the HPS sample, 56 (35%) have [01] /Ha > 0.25 and 106 (65%) 
have [OI]/Ha < 0.25. The corresponding fractions for our sample are nearly identical: 17 
and 34 of our 51 LLAGN fit our definitions strong and weak- [01] nuclei respectively. 

Pig. 1. — Distances and morphological T types for the strong- [01] (triangles, [OI]/Ha > 0.25) 
and weak- [01] LLAGN (circles, [01] /Ha < 0.25) in the HPS sample. Pilled symbols and filled 
areas in the histograms indicate objects in our sample. (All data extracted from HPS97.) 
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Table 1. Sample Properties 



Galaxy 


Spectral 


Hubble 


T-type 


V 


dist. 


pc/ 


log LjJq 


[OI]/Ha 


Name 


Class 


Type 


km 8 


Mpc 


pc 






NGC 266 


LI. 9 SB(rs)ab 


2.0 


4662 


62.4 


302 


39.30 


0.280c 


NGC 315 


LI. 9 E + ; 


-4.0 


4935 


65.8 


319 


39.55 


0.590 


NGC 404 


L2 SA(s)0-: 


-3.0 


-46 


2.4 


11 


37.63 


0.170 


NGC 410 


T2 


E + : 


-4.0 


5295 


70.6 


342 


39.43 


().()97u 


NGC 428 


L2/T2 


SAB(s)m 


9.0 


1160 


14.9 


72 


36.98L 


0.190U 


NGC 521 


T2/H 


SB(r)bc 


4.0 


5039 


67.0 


324 


39.16 


0.086c 


NGC 660 


T2/H 


SB(s)a pec 


1.0 


852 


11.8 


57 


38.89 


0.047 


NGC 718 


L2 SAB(s)a 


1.0 


1734 


21.4 


103 


38.45 


0.210 


NGC 772^ 


H/T2 


SA{s)b 


3.0 


2468 


32.6 


158 


39.24 


().()26b 


NGC 841 


LI. 9 


(R')SAB(s)ab 


2.1 


4639 


59.5 


288 


39.11L 


().5S()b 


NGC 1052 


LI. 9 E4 


-5.0 


1499 


17.8 


86 


39.45r 


0.710 


NGC 1161 


T1.9 


SAO 


-2.0 


1940 


25.9 


126 


38.70 


0.140U 


NGC 1169 


L2 SAB(r)b 


3.0 


2386 


33.7 


163 


38.67 


0.320 


NGC 2681 


Ll.£ 


(R')SAB(rs)0/a 


0.0 


688 


13.3 


64 


38.83L 


0.190b 


NGC 2685^ 


S2/T2 


{R)SBO+pcc 


-1.0 


820 


16.2 


78 


38.66 


0.130b 


NGC 2911 


L2 SA(s)0: pec 


-2.0 


3183 


42.2 


204 


39.38r 


0.310 


NGC 3166 


L2 SAB{rs)0/a 


0.0 


1344 


22.0 


106 


39.10 


0.270 


NGC 3169 


L2 SA(s)a pec 


1.0 


1234 


19.7 


95 


39.02 


0.280 


NGC 3226 


Ll.S 


E2: pec 


-5.0 


1321 


23.4 


113 


38.93 


0.590 


NGC 3245 


T2 


SA(r)0? 


-2.0 


1348 


22.2 


107 


39.69 


0.086c 


NGC 3627 


T2/S2 SAB(s)b 


3.0 


727 


6.6 


31 


38.60 


0.130 


NGC 3705 


T2 SAB(r)ab 


2.0 


1018 


17.0 


82 


38.66 


0.079b 


NGC 4150 


T2 SA(r)0? 


-2.0 


43 


9.7 


47 


38.18 


0.130 


NGC 4192 


T2 SAB(s)ab 


2.0 


-142 


16.8 


81 


38.97 


0.140 


NGC 4438 


LI. 9 SA(s)0/a: 


0.0 


64 


16.8 


81 


39.37 


0.270 


NGC 4569 


T2 SAB(rs)ab 


2.0 


-235 


16.8 


81 


40.28r 


0.062 


NGC 4736 


L2 {R)SA(r)ab 


2.0 


307 


4.3 


20 


37.81r 


0.240 


NGC 4826 


T2 (R)SA(rs)ab 


2.0 


411 


4.1 


19 


38.87r 


0.073 


NGC 5005 


Ll.£ 


SAB(rs)bc 


4.0 


948 


21.3 


103 


39.47 


0.650 


NGC 5055 


T2 SA(rs)bc 


4.0 


504 


7.2 


34 


38.62r 


0.170U 


NGC 5377 


L2 {R)SB{s)a 


1.0 


1792 


31.0 


150 


39.18 


0.250 


NGC 5678 


T2 SAB{rs)b 


3.0 


1924 


35.6 


172 


39.19 


0.079 


NGC 5879 


T2/L2 SA{rs)bc? 


4.0 


772 


16.8 


81 


38.32 


0.160 


NGC 5921 


T2 SB(r)bc 


4.0 


1479 


25.2 


122 


39.16r 


0.110 


NGC 5970 


L2/T2 


SB(r)c 


5.0 


1966 


31.6 


163 


38.06 


0.180c 


NGC 5982 


L2: 


E3 


-5.0 


2904 


38.7 


187 


38.46c 


0.490U 


NGC 5985 


L2 SAB{r)b 


3.0 


2518 


39.2 


190 


38.94 


0.300 


NGC 6340 


L2 SA(s)0/a 


0.0 


1207 


22.0 


106 


38.50 


0.43()b 


NGC 6384 


T2 SAB(r)bc 


4.0 


1667 


26.6 


128 


38.12 


().15()u 


NGC 6482 


T2/S2: 


B: 


-5.0 


3921 


52.3 


263 


39.23 


O.lSOu 


NGC 6500 


L2 SAab: 


1.7 


2999 


39.7 


192 


40.31 


0.230 


NGC 6501 


L2: 


SAO+: 


-0.5 


2869 


39.6 


191 


38.01c 


0.870U 


NGC 6503 


T2/S2 


SA(s)cd 


6.0 


42 


6.1 


29 


37.56 


0.080 


NGC 6702 


L2: 


B: 


-5.0 


4712 


62.8 


304 


38.61c 


0.620U 


NGC 6703 


L2: 


SAO- 


-2.5 


2364 


35.9 


174 


38.51 


0.360U 


NGC 6951^ 


S2 SAB(rs)bc 


4.0 


1424 


24.1 


116 


39.07 


0.230 


IN 1 L 1 1 


T2 SAB{r)b 


3.0 


1147 


18.2 


88 


38.94 


0.140 


NGC 7217 


L2 (R)SA(r)ab 


2.0 


945 


16.0 


77 


38.86 


0.250 


NGC 7331 


T2 SA(s)b 


3.0 


821 


14.3 


69 


38.49 


0.097U 


NGC 7626 


L2: 


E: pec 


-5.0 


3422 


45.6 


221 


38.81 


0.220U 


NGC 7742 


T2/L2 SA{r)b 


3.0 


1653 


22.2 


107 


39.07 


0.130 


NGC 3367 


H SB{rs)c 


5.0 


3041 


43.6 


211 


40.98 


0.031 


NGC 6217 


H {R)SB{rs)bc 


4.0 


1362 


23.9 


115 


40.41 


0.032 


NGC 205 


normal dE5 pec 


-5.0 


-239 


0.7 


3 


< 34.79r 




NGC 221 


normal E2 


-6.0 


-205 


0.7 


3 


< 36.16r 




NGC 224 


normal SA(s)b 


3.0 


-298 


0.7 


3 






NGC 628 


normal SA(s)c 


5.0 


655 


9.7 


3 


< 36.69r 




NGC 1023 


normal SB(rs)0- 


-3.0 


632 


10.5 


51 


< 37.82r 




NGC 2950 


normal (R)SB(r)0 


-2.0 


1337 


23.3 


204 


< 38.42r 




NGC 6654 


normal (R')SB(s)0/a 


0.0 


1823 


29.5 


123 


< 38.11r 




Note. — Col 


(1): Galaxy name; Col. (2): 


Spectral class. 


Col. (3): 


Hubble type. Col. (4): Numerical Hubble 


type. Col. (5): 


Radial velocity. Col. (6): Distance. Col. (7) 


Angular scale. 


Col. (8): 


Ha luminosity. Col. (9): 



[OI]/Ha flux ratio. All quantities were extracted from HFS97. Revised values of log Lnai taken from Ho etal. 
(2003), are indicated by an "r" in column 8. See HFS97 for the meaning of other notes on columns 8 and 9. 

"'^'^ Grouped among TOs in this paper. 

^Grouped with LINERs in this paper (see Perez etal. 2000). 
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2.2. Compcirison with the HFS sample 

Although we have drawn objects exclusively from the HFS sample, it is important 
to evaluate whether we have introduced any bias in our sample selection. In Figure 1 we 
compare the distributions of distances (d) and morphological types (T) in these two samples. 
Triangles and circles correspond to strong and weak- [01] LLAGN respectively, and filled 
symbols mark objects which we have observed. The visual impression from the bottom left 
plot is that our galaxies are well mixed within the HFS sample. This is confirmed by the 
histograms projected along both axis of Figure 1, where we see that our galaxies (filled areas) 
follow reasonably well the distributions of d and T in the full HFS sample (solid lines). In 
both samples the distribution of morphological types of weak- [01] nuclei is skewed towards 
somewhat later type hosts than for strong-[OI] nuclei, but with a substantial overlap, as 
discussed by Ho etal. (2003). 

In Figure 2 we compare the G band equivalent width and the [01] /Ha ratio, both 
extracted from HFS97. Again, one sees that our galaxies are well mixed within the HFS 
galaxies. This is confirmed by a statistical analysis. The mean value and rms dispersion of 
W{G band) is 4.6 ± 1.2 for aU LLAGN in HFS97 and 4.1 ± 1.3 in our sample. No biases 
were detected in the [OI]/Hq; ratio either (0.27 ± 0.21 compared to 0.25 ± 0.19 for the HFS 
and our samples respectively) , and similar results apply to other stellar indices and emission 
line ratios tabulated by HFS97. We thus conclude that our sample is representative of the 
full HFS sample in terms of emission line and stellar population properties. 

We note in passing that Figure 2 already gives away an important result of this in- 
vestigation. Most objects with a weak G band are also weak [01] emitters. The "inverted 
L" shape traced by the galaxies in this plot is only spoiled by a streak of 3 or 4 LINERs 
with [OI]/Hq; ~ 0.6 and W(G band) < 4 A. Apart from these objects, the top-left region 
of Figure 2 is remarkably empty. Of the two objects which intrude most in this "zone of 
avoidance", NGC 5195 and IC 239, the latter has only an upper limit on [01] and both have 
uncertain classifications (HFS97). Similar uncertainties affect several other points in Figure 
2, but these do not affect the global shape of the distribution. Other versions of this plot 
are presented and discussed in §6. 

Fig. 2. — Equivalent width of the G band versus [OI]/H« for LLAGN in the HFS sample. 
Filled symbols mark objects in our sample. The horizontal line at [01] /Ha = 0.25 indicates 
our border-line between weak and strong-[OI] nuclei. (All data extracted from HFS97.) 
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3. Observations 

Observations were carried out in four runs at the 2.5m Nordic Optical Telescope and one 
run at the Kitt Peak National Observatory 2.1m. The NOT observations used the ALFOSC 
detector with grism number 6 and a 1" wide slit, which gives a dispersion of 1.4 A per pixel 
and covers the wavelength range 3500-5500 A. The sht was oriented along the paralactic 
angle. Most observations were made under sub-arcsecond seeing and photometric conditions. 
The KPNO observations were carried out using the Gold Camera with grating 26 old and a 
slit 2" wide, which gives a dispersion of 1.24 A per pixel, and covers the 3400-5500 A interval. 
The slit was oriented in the N-S direction and the observations were done as close to the 
meridian passage as possible to avoid light losses by differential refraction. The weather 
during the observations was not photometric, and the seeing was typically 2-3". A log of 
the observations is given in Table 2. 

The data were reduced using standard IRAF procedures. The individual frames were 
overscanned, bias subtracted and flat fleld divided. The calibration also followed standard 
steps. The wavelength scale was calibrated using an HeNeAr lamp. The flux calibration 
was done with observations of spectrophotometric standard stars observed with a slit width 
of 10". The individual frames were combined using crreject task in IRAF to remove cosmic 
rays. Sky was subtracted with the background task fitting a polynomium along the spatial 
hght profile. 

3.1. Nuclecir Extractions 

Most of the spectra have excellent signal out to several arcsec from the nucleus. In Cid 
Fernandes etal. (2003b, Paper III) we explore this spatial information to map the stellar 
populations in the central regions of these objects, in analogy with previous work for type 
2 Seyferts (Cid Fernandes, Storchi-Bergmann & Schmitt 1998; Joguet etal. 2001; Raimann 
etal. 2003). In this paper we concentrate on the nuclear spectra. Off-nuclear extractions 
will only be used as templates for the nuclear spectra (§4). 

Nuclear spectra were extracted from the central l"xl.l" (6 pixels) for the 52 galaxies 
observed at NOT and 2"x2.3" (3 pixels) for the 8 galaxies observed at KPNO. Using the 
distances in Table 1, the projected areas covered by these extractions correspond to radii of 
11-204 pc, with a median value of 71 pc considering only the LLAGN in the sample. The 
signal-to-noise ratio was estimated from the rms fluctuation in the 4789-4839 A interval, 
which is free of major absorption or emission features. We obtain S/N in the 23-88 range, 
with a median value of 50. In the 4010-4060 A interval the median S/N is 22. 
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Table 2. Log of observations 



Name 


Telescope 


Date 


Exposure (s) 


airmass 


NGC 266 


NOT 


8-21-01 


3x1200 


1.03 


NGC 315 


NOT 


8-22-01 


3x1200 


1.11 


NGC 404 


KPNO 


7-26-01 


2x1200 


1.29 


NGC 410 


NOT 


8-21-01 


3x1200 


1.13 


NGC 428 


NOT 


8-23-01 


4x1200 


1.26 


NGC 521 


NOT 


8-23-01 


3x1200 


1.69 


NGC 660 


NOT 


8-21-01 


4x1200 


1.07 


NGC 718 


NOT 


8-22-01 


3x1200 


1.20 


NGC 772 


NOT 


8-23-01 


3x1200 


1.04 


NGC 841 


NOT 


8-22-01 


4x1200 


1.02 


NGC 1052 


NOT 


8-22-01 


2x1200 


1.26 


NGC 1161 


NOT 


8-23-01 


2x1200 


1.06 


NGC 1169 


NOT 


11-8-02 


3x1200 


1.05 


NGC 2681 


NOT 


11-8-02 


3x900 


1.33 


NGC 2685 


NOT 


11-8-02 


3x1200 


1.03 


NGC 2911 


NOT 


5-3-03 


3x1200 


1.12 


NGC 3166 


NOT 


6-4-03 


2x1200 


1.11 


NGC 3169 


NOT 


5-3-03 


3x1200 


1.14 


NGC 3226 


NOT 


5-14-01 


3x1200 


1.05 


NGC 3245 


NOT 


5-12-01 


1800-1-1200+600 


1.04 


NGC 3627 


NOT 


5-13-01 


3x1200 


1.09 


NGC 3705 


NOT 


5-14-01 


3x1200 


1.12 


NGC 4150 


NOT 


5-14-01 


4x1200 


1.28 


NGC 4192 


NOT 


5-12-01 


4x900 


1.04 


NGC 4438 


NOT 


5-15-01 


4x1200 


1.25 


NGC 4736 


NOT 


5-3-03 


3x900 


1.04 


NGC 4569 


NOT 


5-12-01 


3x1200 


1.11 


NGC 4826 


NOT 


5-14-01 


3x1200 


1.26 


NGC 5005 


NOT 


5-3-03 


3x1200 


1.09 


NGC 5055 


NOT 


5-16-01 


4x1200 


1.16 


NGC 5377 


NOT 


6-3-03 


3x1200 


1.24 


NGC 5678 


NOT 


5-13-01 


4x1200 


1.18 


NGC 5879 


NOT 


5-14-01 


4x1200 


1.21 


NGC 5921 


NOT 


5-12-01 


3x1200 


1.09 


NGC 5970 


NOT 


8-23-01 


3x1200 


1.19 


NGC 5982 


NOT 


6-4-03 


3x1200 


1.28 


NGC 5985 


NOT 


8-22-01 


3x1200 


1.29 


NGC 6340 


KPNO 


7-26-01 


3x1200 


1.34 


NGC 6384 


NOT 


5-13-01 


3x1200 


1.12 


NGC 6482 


NOT 


5-12-01 


3x1200 


1.03 


NGC 6500 


NOT 


5-12-01 


3xl200-f900 


1.02 


NGC 6501 


NOT 


8-22-01 


3x1200 


1.10 


NGC 6503 


NOT 


5-13-01 


3x1200 


1.33 


NGC 6702 


NOT 


8-21-01 


3x1200 


1.06 


NGC 6703 


KPNO 


7-24-01 


3x1200 


1.03 


NGC 6951 


KPNO 


7-24-01 


2x1200 


1.21 


NGC 7177 


KPNO 


7-25-01 


3x1200 


1.20 


NGC 7217 


KPNO 


7-26-01 


3x1200 


1.01 


NGC 7331 


KPNO 


7-24-01 


3x1200 


1.03 


NGC 7626 


NOT 


8-23-01 


3x1200 


1.40 


NGC 7742 


KPNO 


7-24-01 


2x1200 


1.10 


NGC 3367 


NOT 


5-13-01 


3x1200 


1.04 


NGC 6217 


NOT 


8-23-01 


3x1200 


1.62 


NGC 205 


NOT 


11-7-02 


3x900 


1.03 


NGC 221 


NOT 


11-7-02 


2x1200 


1.22 


NGC 224 


NOT 


11-7-02 


2x900-H600 


1.08 


NGC 628 


NOT 


11-7-02 


3x1200 


1.03 


NGC 1023 


NOT 


11-7-02 


3x900 


1.10 


NGC 2950 


NOT 


11-7-02 


3x900 


1.21 


NGC 6654 


NOT 


5-16-01 


3x1200 


1.41 
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3.2. HFS extractions 

Most of the HFS97 spectra were collected through a 2" x 4" aperture, ~ 7 times larger 
in area than our typical nuclear extractions. In order to allow a more meaningful compar- 
ison of our measurements with those in HFS97 we have also extracted spectra through a 
2" X 4" aperture (hereafter the "HFS-aperture"). This was straightforward for the KPNO 
observations, which have a slit width equal to that in HFS97. In order to construct an 
HFS-aperture for the NOT data, which was collected through a 1" slit, we first extracted 
a 4" long spectrum and then multiplied it by 2 except in the nucleus (ie, the inner 1.1"), 
which was represented by an interpolation between off- nuclear extractions centered at ±1.1" 
to avoid overweighting any nuclear point source. The flux in this approximate HFS-aperture 
is typically 4 times larger than in our nuclear extractions. 

All spectral properties discussed in §5 were measured both through our nuclear extrac- 
tions and through this HFS-aperture, but only the nuclear values are published here, since 
except for galaxies with noticeable gradients (discussed in Paper III) the results are similar 
for both apertures. Tables for HFS-apertures measurements are available upon request. 

Besides spatial resolution, our observations differ from those of HFS97 in spectral cov- 
erage. The blue spectra in HFS97 cover the 4230-5110 A interval, while our spectra go from 
3500 up to 5500 A. As weU shaU soon see, the region between 3500 and 4200 A proved to be 
most revealing. 

4. Spectra 

In this section we present the spectra and two different methods to characterize the stel- 
lar populations in the central regions of LLAGN. The first method consists of an empirical, 
qualitative classification based on the comparison with non-active galaxies representative of 
simple stellar populations with different ages. The second method quantifies this classifica- 
tion by modeling the starlight using the non-active galaxies as spectral templates. In Section 
5 we present measurements of several spectral indices (such as the 4000 A break, R5, Call 
K, G band and others), which provide yet another way to characterize stellar populations. A 
comparison of all these methods is important in view of the diversity of methods to analyse 
the stellar populations of AGN found in the literature. This comparison allows us to do a 
relative calibration between them, and can be used to convert the results obtained with one 
method to the others. 



-13- 



4.1. Normal galctxies and an empirical stelleir population classification 

We start our presentation of the nuclear spectra by the 9 non-active galaxies in our 
sample, shown in Figure 3. The order of the spectra in this plot reflects approximately 
the age of the dominant stellar population, younger systems at the top. These normal 
galaxy spectra are used in this work as a reference basis for an empirical stellar population 
classification. Galaxies in Figure 3 can be separated into four classes: 

(y) Galaxies with evidently young, < 10^ yr, stellar populations (NGC 3367 and NGC 
6217), characterized by a blue continuum and weak metal absorption lines; 

(/) galaxies with a dominant intermediate age, 10^-10^ yr population, characterized by 
pronounced HOBLs (NGC 205); 

{I/O) galaxies with a mixture of intermediate age and older stars (NGC 221 and NGC 628); 

(O) galaxies dominated by an old stellar population (NGC 224, NGC 1023, NGC 2950 and 
NGC 6654), characterized by strong metal hnes. 

For short, these four stellar population classes are denominated Y, 7, I/O and O, 
respectively. We use the symbol 77 to denote these classes. 

We have compared each of the LLAGN in our sample to these normal galaxies and 
classified them onto these same 4 stellar population categories according to a similarity 
criterion. This was first done by a simple visual inspection of the spectra, and later confirmed 
by means of the starhght modehng scheme described in §4.4. The results are listed in the 
last column of Table 3. None of the objects resembles predominantly young systems like 
NGC 3367, so all our LLAGN fall onto the 77 = /, //O and O categories. 

This empirical classification provides a stellar population counterpart to the LINER/TO 
classification, which is based on emission hne properties. The 77 classes are used here just as 
a heuristic aid to sort the LLAGN spectra according to their main stellar populations. This 
is a useful strategy to investigate connections between emission line and stellar population 
properties. However, classification criteria, while certainly helpful to identify general trends, 
have a certain degree of uncertainty and arbitrariety. When discussing our results we will thus 

Fig. 3. — Spectra of the normal galaxies in the sample, including two Starburst nuclei (top). 
Galaxies are sorted according to the age of the dominant stellar population, with younger 
systems at the top. F, /, I/O and O denote our empirical stellar population classes. 
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not restrict ourselves to comparisons between discrete emission line and stellar population 
classes. In §6 we give more emphasis to comparisons between continuous stellar population 
and emission line properties. We shall see that a clearer picture emerges if one avoids the 
rigidity of taxonomy. 

4.2. LINERs and Transition Objects 

The nuclear spectra for all LLAGN in our sample are shown in Figures 4-8. Since our 

main aim in this paper in to investigate the relation between the stellar populations of these 
objects and their emission line properties, we have split the spectra in six slots: strong- [01] 
sources with either O, I/O and / stellar population characteristics (Figures 4 and 5) and 
weak-[OI] with either rj = O, I/O and / (Figures 6, 7 and 8). 

As is typical of LLAGN, the spectra are completely dominated by starlight. With few 
exceptions (eg, NGC 1052 and NGC 6501), emission lines are generally weak, and sometimes 
altogether absent (eg, NGC 410 and NGC 5055). Important diagnostic lines like H/? and 
[0111] AA4959, 5007 arc rarely seen, although [OII]A3727 is detected in most spectra. As 
thoroughly discussed and illustrated by HFS97, starhght subtraction is an absolute necessity 
in order to measure emission lines in these objects, so we postpone the analysis of emission 
lines to a future communication. Except for [Oil], lines fluxes for galaxies in the present 
sample have already been measured by HFS97, including the important region around Ha, 
so we will use their measurements when necessary. In the following we concentrate on the 
measurement and analysis of stellar population features. 

As is evident from Figures 4-8, there is a clear connection between stellar population 
characteristics and the emission line properties encoded in the weak/strong- [01] classifica- 
tion. Old stellar populations are substantially more frequent in strong than in weak- [01] 
sources, whereas the opposite happens for stellar population class /. Numerically, we find 
that out of our 17 strong- [01] nuclei, 10 (59%) are in the rj = O class, 5 (29%) have rj = I/O 
and only 2 (12%) belong to the rj = I class, while the corresponding numbers for the 34 
weak-[OI]'s are 7 (21%), 11 (32%) and 16 (47%) respectively. 

Fig. 4. — Nuclear spectra of strong- [01] LLAGN with predominantly old stellar populations 
{V = 0). 



- 15 - 



4.3. High Order Balmer Absorption Lines 

The most notable feature in this data set is the ubiquity of the HI high order Balmer 
series in absorption (H8 and higher), particularly among weak- [01] nuclei. By definition, 
objects where these features are clearly present are attributed rj = I (ie, they resemble NGC 
205). Systems with HOBLs are also bluer and have shallower metal absorption lines than 
the rest. Overall, their spectra reveal clear signs of a 10^-10^ yr intermediate age, "post- 
starburst" population. HOBLs may also be present in I/O nuclei, but blended with metal 
lines of old stars to such an extent that they do not stand out clearly in the total spectrum. 

HOBLs have been seen before in a few LLAGN, the most famous example being NGC 
4569 (Keel 1996), which is also in our sample (Figure 8a). However, an incidence rate 
as high as that suggested by our survey had not been suspected before. Ho etal. (2003) 
argue that objects like NGC 4569 are not representative of LLAGN in their sample, most 
of which have a "demonstrably old" stellar population. While we agree that NGC 4569 
is an extreme example and that old stars dominate the light in most LLAGN, our spectra 
reveal the presence of intermediate age stars in at least 1/3 of the objects. The reason why 
we identify more of these populations than Ho et al. (2003) is that they had to base their 
judgment on the strength of the H7, H/9 and Ha absorptions, which are much more affected 
by emission and dilution by old stars than HOBLs. 

As remarkable as the high frequency of HOBLs in weak- [01] 's is their apparent dearth 
among strong- [01] nuclei. HOBLs are ~ 5 times more frequent in weak than in strong- [01] 
sources [r] = I fractions of ~ 50% compared to ~ 10% respectively). These findings reveal a 
strong link between the central stellar populations and emission line properties, particularly 
[Olj/Ha. 

The statistics above obviously depend on our adopted borderline between weak and 
strong-[OI] nuclei (at [OI]/Hq; = 0.25), but the dichotomy persists (albeit somewhat diluted) 
even if one adopts the HFS97 [Olj/Ho; = 0.17 classification criterion. A less classification 
dependent assessment of the relation between stellar and emission line properties is presented 
in §6, where we replace the discrete stellar population and emission line classes used above 
by a more continuous description. 

Fig. 5. — Nuclear spectra of strong- [01] LLAGN with stellar population classes rj — I/O (a) 
and rj — I (b) . 
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4.4. Empirical steirlight modeling 

As noted in §1, the location of TOs in between LINERs and Hll-nuclei in diagnostic 
diagrams suggests that their emission hne spectrum results from a mixture of AGN spectrum 
plus photoionization by young stars (HFS97). If this is the case, one would expect to find 
direct signatures of stellar populations of 10^ yr or less in their spectra. Although the high 
incidence rate of HOBLs among weak- [01] nuclei points to some sort of connection between 
stellar populations and emission line properties, these features signal the presence of a ~ 10^- 
10^ yr post-starburst-like component, not a young starburst. The clearest single signature 
of young stars in the spectral range covered by our data is the WR bump at 4680 A. This 
feature is not evident in any of our LLAGN spectra. The WR bump was only detected in 
the Starburst galaxy NGC 3367 (Figure 3), which is part of our comparison sample. 

Given the dominance of old and intermediate age stars in LLAGN, a weak young star- 
burst (with or without a correspondingly weak WR bump) could well be lurking in some of 
our spectra. Detailed modeling of the stellar population mixture is required to evaluate this 
possibility. This modeling is performed in a separate communication (Paper II), in order not 
to mix data with models. In this paper we want to keep the analysis at a more empirical 
level. We have thus implemented an entirely self contained starlight modeling technique 
which, instead of resorting to theoretical spectra, uses only our observed spectra. 

Each spectrum was modeled as a combination of a base of template spectra. Two bases 
were considered: A "normal galaxy base", containing five spectra from our comparison 
sample (NGC 3367, NGC 205, NGC 221, NGC 1023 and NGC 2950), and an "off-nuclear 
base", containing two off-nuclear extractions centered at ±2.4" from the nucleus. Except 
for the possibility of using off-nuclear templates, this procedure is very similar to the one 
employed by HFS97 in their starlight subtraction scheme. 

All spectra were normalized to the flux at 4020 A and corrected for Galactic extinction 
prior to the modeling. The Ab values of Schlegel, Finkbeiner Davis (1998), extracted from 
NED^, and the reddening law of CardeUi, Clayton & Mathis (1989) were used for this purpose. 
Intrinsic extinction by an uniform screen of dust was allowed for. An algorithm was developed 
which searches for the combination of base spectra and extinction which best reproduces the 



^Thc NASA/IPAC Extragalactic Database (NED) is operated by the .Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 

Fig. 6. — Nuclear spectra of weak-[OI] LLAGN with stellar population classes 77 = O (a) 
and 77 = I/O (b). 
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nuclear spectrum of each LLAGN. The code was adapted from the empirical population 
synthesis code of Cid Fernandes et al. (2001b) by simply replacing its spectral base by the 
full spectra in either one of these two bases. Regions around emission lines were masked 
out in the comparison between observed and model spectra. For the off-nuclear base, in 
most cases it was sufficient to mask out regions of 50 A around [Oil], H/3 and [0111]. For 
the normal galaxy base we further masked regions around H7, B.6 and [N1]A5200 since these 
are rather strong in the spectrum of NGC 3367, included in this base to represent a young 
star burst. 

The presence of an emission line object in the normal galaxy base is another differ- 
ence with respect to the starlight modeling scheme employed by HFS97. This difference 
is easily understood. HFS97 were interested in measuring emission lines from a starlight 
free spectrum, so their starlight templates were constructed out of galaxies without emission 
lines. Our interest here is quite the opposite: We want to analyze the starlight itself, not 
the emission lines. NGC 3367 was introduced in the fitting to evaluate how strong a young 
starburst can be accommodated in our LLAGN spectra. A consequence of this choice is that 
our observed minus model residual spectra will have underestimated, emission line fluxes 
when the NGC 3367 component is present in a significant proportion. There is no way of 
circumventing this without editing out the emission lines from NGC 3367 or resorting to 
theoretical stellar population spectra, since an emission line free starburst does not exist 
in nature. Nevertheless, in the majority of cases we find that the NGC 3367 component 
contributes very httle (< 2%) to the optical flux. In these cases, our residual spectra can be 
seen as pure emission spectra. 

Examples of the results obtained with both spectral bases are shown in Figure 9. The 
top panels show the observed spectrum (F^^^, thick hne) and the two models for the starlight 
spectrum. The middle spectrum shows F^'^^^ the template constructed out of linear combi- 
nations of the normal galaxy base, while the bottom spectrum is the combination 
of two off-nuclear extractions which best matches the nuclear spectrum. Both models were 
shifted vertically for clarity. The residual spectra are shown in the bottom plots. Panels 
d-f correspond to F^^^ — p^^"^^P^ while panels g-i show F^^^ — Except for regions 
containing emission lines, the typical difference between F^'"^ and the i^^'^^P model is better 
than 3% for all cases shown. This is comparable to the noise level. Off-nuclear templates 
produce equally good results for the examples shown, with residuals of 4% for NGC 266 and 
NGC 6500 and 1% for NGC 2681. In some cases (eg, NGC 772, NGC 4569, NGC 5921) 

Fig. 7. — Nuclear spectra of weak- [01] LLAGN with stellar population class rj = I/O (a) 
and 77 = / (b). 
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the F^^^ — residual reveals a nuclear blue component not present in the off-nuclear 

spectra. This happens specially in objects which have strong young and/or intermediate age 
populations as revealed by the strength of the NGC 3367 and NGC 205 components in the 
normal galaxy decomposition. We also note in passing that Figure 9 shows that HOBLs arc 
present not only in the nucleus of NGC 2681, but also in off-nuclear extractions. The region 
producing the HOBLs is therefore spatially extended, a result which is examined in more 
detail in Paper III. 

One thus sees that, except for galaxies with noticeable gradients, both methods yield 
comparable accuracy. For our current purposes, the normal galaxy base is more relevant, 
since it provides a way of quantifying the stellar population mix in LLAGN. In fact, this de- 
composition method provides a quantitative basis for the comparative stellar population clas- 
sification outlined in §4.1. To illustrate how this works, one can write the relative contribu- 
tions at A = 4020 A of NGC 3367 (rj^Y), NGC 205 (r? = /), NGC 221 {rj = I/O) and NGC 
2950 plus NGC 1023 {rj = O) as a normalized four-component vector x = {xy,xj,Xi/o,Xo)- 
For the ri = nuclei of NGC 266 and NGC 7626 we find x = (0, 5, 18, 77)% and (1, 2, 0, 97)% 
respectively, while for NGC 4438, which we classified ast] = I/O, this vector is (0, 3, 51, 46)%, 
and for the rj = I nucleus of NGC 5921 we find (2,70,4,24)%. In over 90% of the cases 
this scheme yields the same classification as that initially inferred by visual inspection of the 
spectra. In the few cases where a disagreement occurred, we adopted the visually estimated 
77 class, although this has no consequence for the results of this paper. 

The results of the starlight modeling in terms of the normal galaxy base are summarized 
in Table 3. A first important result of this analysis is that 23 of our 51 LLAGN contain a 
NGC 205 component stronger than 20%. Except for NGC 5879, which has a noisy spectrum, 
all nuclei where xj > 30% were classified as rj — I because of their conspicuous HOBLs. 
Between xj — 10% and 30% some objects were classified as rj — I (NGC 772) and others as 
I/O (eg, NGC 4192, NGC 5055). The starlight analysis therefore confirms the high incidence 
rate of HOBLs in LLAGN, specially weak-[OI]'s, and helps quantifying the presence of these 
intermediate age populations in cases where their contribution to the optical spectrum is not 
visually obvious. 

Fig. 8. — Nuclear spectra of weak- [01] LLAGN with a strong intermediate age stellar pop- 
ulation (r) — I). 
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Table 3. Starlight modeling results 



Galaxy 


Xy 


XI 




xo 


Av 


A 






(NGG 3367) 


(NGC 205) 


(NGG 221) 


(NGG 1023 + NGC 2950) 








NGC 0266 


0.1 


5.2 


17.8 


77.0 


0.18 


2.8 


O 


NOG 0315 


3.8 


9.2 


1.5 


85.5 


0.09 


3.5 


O 


NGG 0404 


0.0 


78.5 


0.0 


21.5 


0.90 


6.6 


I 


NGG 0410 


0.0 


1.8 


0.7 


97.5 


0.04 


2.6 


o 


NGG 0428 


0.0 


49.8 


0.0 


50.2 


0.31 


8.7 


I 


NGC 0521 


1.3 


11.1 


0.1 


87.6 


0.31 


2.7 


o 


NGG 0660 


0.0 


59.4 


0.0 


40.6 


2.38 


11.7 


I 


NGG 0718 


0.0 


51.3 


2.3 


46.4 


0.13 


2.7 


I 


NGG 0772 


32.1 


10.9 


24.4 


32.5 


0.38 


1.8 


I 


NGG 0841 


0.0 


37.7 


17.4 


44.9 


0.18 


2.3 


I 


NGC 1052 


16.5 


0.0 


0.5 


83.0 


0.57 


3.0 


o 


NGC 1161 


0.0 


6.7 


0.1 


93.2 


0.13 


2.7 


o 


NGC 1169 


0.0 


0.0 


48.2 


51.8 


0.26 


5.7 


I/O 


NGG 2681 


0.0 


66.2 


1.4 


32.4 


0.22 


3.4 


I 


NGG 2685 


0.3 


0.9 


50.5 


48.4 


0.21 


4.9 


I/O 


NGC 2911 


0.0 


3.4 


0.0 


96.6 


0.83 


4.6 


o 


NGC 3166 


0.0 


28.8 


9.4 


61.8 


0.19 


2.3 


I/O 


NGC 3169 


0.0 


21.1 


0.0 


78.9 


0.73 


4.4 


I/O 


NGC 3226 


0.0 


2.7 


0.0 


97.3 


1.08 


4.2 


o 


NGG 3245 


14.4 


3.8 


25.5 


56.3 


0.17 


1.5 


I/O 


NGG 3627 


0.0 


66.6 


1.4 


32.0 


0.47 


3.0 


I 


NGG 3705 


0.0 


39.3 


13.8 


46.9 


0.54 


3.6 


I 


NGC 4150 


0.0 


66.8 


32.0 


1.2 


0.00 


6.5 


I 


NGC 4192 


0.0 


19.1 


34.7 


46.1 


1.41 


3.4 


I/O 


NGC 4438 


0.0 


2.6 


50.9 


46.5 


1.48 


2.9 


I/O 


NGG 4569 


28.2 


71.8 


0.0 


0.0 


0.00 


5.3 


I 


NGG 4736 


0.0 


48.8 


0.0 


51.2 


0.07 


2.3 


I 


NGG 4826 


0.0 


32.8 


10.7 


56.5 


0.24 


2.3 


I 


NGG 5005 


0.0 


45.1 


0.4 


54.5 


0.77 


3.7 


I 


NGC 5055 


2.0 


20.8 


9.1 


68.1 


0.24 


3.4 


I/O 


NGC 5377 


0.0 


76.8 


0.1 


23.1 


0.02 


2.5 


I 


NGC 5678 


0.0 


89.7 


0.3 


9.9 


1.22 


3.8 


I 


NGG 5879 


0.0 


52.0 


0.9 


47.1 


0.46 


7.3 


I/O 


NGG 5921 


1.7 


70.1 


3.7 


24.5 


0.00 


2.8 


I 


NGG 5970 


0.0 


5.7 


56.4 


37.9 


0.13 


4.6 


I/O 


NGC 5982 


0.2 


10.3 


2.5 


87.0 


0.02 


1.9 


o 


NGC 5985 


0.0 


6.8 


34.4 


58.8 


0.19 


4.9 


o 


NGC 6340 


0.0 


4.1 


0.0 


95.9 


0.98 


6.6 


o 


NGC 6384 


1.2 


2.2 


61.3 


35.3 


0.02 


4.4 


I/O 


NGG 6482 


10.1 


0.2 


0.0 


89.7 


0.12 


2.5 


o 


NGG 6500 


24.7 


0.1 


4.0 


71.3 


0.08 


2.3 


I/O 


NGC 6501 


1.2 


11.5 


0.0 


87.3 


0.25 


2.3 


o 


NGC 6503 


22.5 


55.8 


16.3 


5.4 


0.32 


3.0 


I 


NGC 6702 


0.0 


7.0 


18.6 


74.4 


0.00 


2.7 


o 


NGG 6703 


0.0 


3.4 


0.0 


96.6 


0.22 


3.5 


o 


NGG 6951 


2.0 


29.4 


18.9 


49.7 


0.41 


3.1 


I/O 


NGG 7177 


2.4 


13.7 


43.1 


40.8 


0.39 


2.8 


I/O 


NGG 7217 


0.0 


7.2 


0.0 


92.8 


0.43 


3.6 


o 


NGC 7331 


0.0 


5.3 


15.0 


79.8 


0.37 


2.8 


o 


NGC 7626 


1.1 


1.8 


0.0 


97.2 


0.09 


2.2 


o 


NGC V712 


0.1) 


11.:-! 


D.O 




O.K) 


1.1 


I/O 



Note. — Results of the starlight decomposition in terms of the normal galaxies NGC 3367, NGC 205, NGC 
221, NGG 1023 and NGC 2950. Cols. (2)-(5); xy, xi, Xj/o £^nd xo, given as percentage fractions of the flux 
at 4020 A. Gol. (6): V band extinction (in mag) of the best fit. Col. (7): Mean absolute percentage difference 
between the observed and model spectra: A =< IF^*"^ - Fmodeijy^obs > ^g). stellar population class. 
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4.5. The weEikness of young stellsir populations in LLAGN 

A second result of the starlight modeling is that few LLAGN have an optically relevant 
young starburst. The contribution of the NGC 3367 component is larger than 20% in just 
four cases. In order of decreasing xy, these are NGC 772 (32%), NGC 4569 (29%), NGC 
6500 (25%) and NGC 6503 (22%), which are all weak- [01] 's. 

Of all our LLAGN, only NGC 6500 presents marginal evidence for the presence of the 
WR bump. This can be seen by the coherent broad residual with amplitude of no more than 
a few % in the F^^^ — f^^^^ residual spectrum in Figure 9f. Naturally, the amplitude of this 
residual increases a bit removing the WR bump from the NGC 3367 base component, but 
even then we cannot claim a conclusive detection of this feature. One way to double check 
this possibility is to evaluate the residual spectrum obtained with the off-nuclear base. If 
WR stars are concentrated in the nucleus, a WR bump should appear in the F^^^ — p°^~"^'"^ 
difference spectrum. In the case of NGC 6500, this residual reveals a weak narrow nebular 
HeIIA4686, but no clear sign of a broad bump (Figure 9i). A tentative detection of broad 
Hell in this galaxy has been previously reported by Barth etal. (1997). 

Detecting the WR bump is already hard in bona fide Starburst galaxies (eg, Schaerer, 
Contini & Kunth 1999), and the superposition of a dominant older stellar population only 
makes matters worse. In NGC 3367 this feature has an equivalent width of ~ 5 A, similar 
to the values in WR galaxies (Schaerer, Contini & Pindao 1999) and metal rich giant HII 
regions (Pindao etal. 2002). Hence, a NGC 3367-like component diluted to < 30% of the 
optical flux yields WR bump equivalent widths of at most 1.5 A, which is comparable to 
our spectral resolution, so weaker features are unhkely to be detected. The weakness of the 
young starburst component in LLAGN may thus explain why we have not detected the WR 
bump in our data. 

It is therefore clear that young stellar populations contribute very little to the optical 
spectra of LLAGN. It is important to point out that even when < 10^ yr stars dominate 
the UV light, as in NGC 404, NGC 4569, NGC 5055 and NGC 6500 (Maoz etal. 1998), 



Fig. 9.— Results of the starhght modehng for NGC 266 (77 = O, strong-[OI]), NGC 2681 
(77 = 7, weak[OI]) and NGC 6500 (77 = I/O, weak-[OI]). (a-c): The top panels show the 
nuclear spectrum F^'^^ (top) and two starhght models: F*^™p (middle), constructed out of 
hnear combinations of the normal galaxies NGC 3367, NGC 205, NGC 221, NGC 2950 and 
NGC 1023; and p"^^'^^^'^ (bottom), constructed combining two off- nuclear extractions. All 
spectra are normalized to the flux at 4020 A, and the two models are shifted downwards for 
clarity, (d-f): F"^^ - FI^'^p residual spectra, (g-i): ^ - Ff 
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their contribution does not exceed 30% of the A4020 flux. In fact, these young starbursts 
are not even detected in the optical for NGC 404 and NGC 5055 due to their low contrast 
with respect to older populations. 

While current star formation seems to be proceeding at a residual level in LLAGN, the 
ubiquity of intermediate age populations indicates that its was much more prominent 10^- 
10^ yr ago, which points to a time-decaying star formation activity. This result contrasts 
with that obtained for studies of type 2 Seyferts, where young stellar populations are found 
in abundance (eg, Cid Fernandes etal. 2001a; Joguet etal. 2001). The mean age of stars in 
the central regions of Scyfcrt 2s is thus smaller than that in LLAGN. It is equally instructive 
to compare the stellar populations of weak and strong- [01] LLAGN. The finding that inter- 
mediate age populations are found almost exclusively among weak- [01] nuclei implies that 
these sources are on-average younger than those with strong [01] emission. The implications 
of these results for evolutionary scenarios are considered in §6 and Paper II. 

5. Spectral Properties 

In this section we present measurements of a suite of stellar features in the nuclear 
spectra presented above. Measurements through the HPS aperture discussed in §3.2 are 
only used here for comparison purposes. Tables with these and other measurements not 
published here are available upon request. 

5.1. Equivalent widths and colors in Bica's system 

In a series of papers starting in the mid 80's, Bica & co-workers have explored a stellar 
population synthesis technique known as Empirical Population Synthesis (EPS), which de- 
composes a given galaxy spectrum in a combination of observed spectra of star clusters. In 
practice, instead of modehng the spectrum directly, this method synthesizes a number of 
absorption line equivalent widths (Wx) and continuum colors {C\), which are used as a com- 
pact representation of Fx. EPS has proven a very useful tool in the analysis of galaxy spectra 
in a variety of contexts, from normal galaxies (Bica 1988), to starbursts (Raimann et al. 2000; 
Cid Fernandes, Leao, Rodrigues-Lacerda 2003) and even AGN (Schmitt, Storchi-Bergmann 
& Cid Fernandes 1999; Raimann etal. 2003). 

In this section we present measurements of Wx and Cx in Bica's system. These will 
be used here as a raw measures of stellar population characteristics in our sample galaxies. 
These same data are used in Paper II as input for an EPS analysis. As a side-step to this 
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investigation, we have developed a fully automated scheme to measure W\ and C\ in this 
system. This is discussed next. 

5.1.1. Automated measure of the Pseudo Continuum 

A drawback of Bica's system is its subjectivity. Both equivalent widths and colors are 
measured relative to a pseudo continuum (PCx), defined at pre-selected pivot wavelengths, 
which is traced interactively over the observed spectrum (see Cid Fernandes et al. 1998 for 
an illustrated discussion). Historically, this has limited the use of EPS to a small number 
of researchers initiated in the "art" of drawing this continuum. Furthermore, the lack of an 
objective recipe to define PCx makes this whole system useless in face of the huge spectral 
databases currently available, which must be processed in an entirely automated manner. 
The present data set is itself an example: Counting all nuclear and off-nuclear extractions, 
there are over 700 spectra. 

In order to remedy this situation, we have formulated an a posteriori definition of Bica's 
system. This was done by means of a straightforward computational procedure which mimics 
the placement of the pseudo continuum, which until now has been carried out by hand. 

We started by measuring by hand a representative sub-set of 42 spectra in the current 
data set, all normalized to the respective median fiux in the 4789-4839 A interval. For each 
pivot wavelength A we have then computed median fluxes in windows of different sizes and 
locations. Windows were restricted to be at least 10 A wide and centered within 100 A of 
the corresponding pivot A. The manually measured PCx fluxes were then compared to the 
automatically measured median fluxes (MFx) for all spectra, and a linear relation 

PCx = ax + bxMFx (1) 

was fltted using an ordinary least squares. The optimal window size and location where cho- 
sen to be those which produced the best linear correlation coefficient (ie, smallest residuals). 
Window parameters and the corresponding MFx to PCx conversion coefficients ax and bx 
are listed in Table 4 for six pivot A's used in this work: 3660, 3780, 4020, 4510, 4630 and 
5313 A. 

This recipe works extremely well, with Hncar correlation coefficients above r = 0.97 for 
all A's. Differences between manual and automated PCx measurements in our training set 
of 42 galaxies are of order 3%. This is well within the level of agreement between manual 
measurements of PCx carried out by different people (Cid Fernandes etal. 1998). We have 
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further verified that the inclusion of other terms in equation (1) does not improve the fits 
significantly. Figure 10 illustrates the pseudo continuum, as measured through our objective 
system. Equivalent widths are measured with respect to the continuum formed by linear 
interpolation between the PC\ pivot points. The differences between W\ measured with 
manual and automatic PCx is also of order 3%, being always smaller than 0.5 A. 

We thus conclude that we have succeeded in formulating a non-subjective definition of 
the pseudo continuum (and hence equivalent widths) in Bica's system, thus overcoming a 
long held objection to this measurement system. 



5.1.2. Results 

We have measured 7 equivalent widths in Bica's system: Wc, Wy,ib, Wk, Wh, Wcn, 
Wq and Wmq- The window definitions for these absorption features arc the same defined in 
previous works (Bica & Alloin 1986; Bica, Alloin & Schmidt 1994). Note that Wc is actually 
centered in a continuum region just to the blue of H9. In old populations, the spectrum in 
this region is well bellow the pseudo continuum due to a multitude of weak metal lines, 
yielding Wc ~ 4-5 A (eg., NGC 266 and NGC 521 in Figure 10). For populations of ~ 1 
Gyr or less, the continuum raises and Wc approaches 0, as seen in NGC 4150 and NGC 4569 
(Figure 10). For our sample, this index is a useful tracer of the presence or HOBLs. All 
cases where the HOBLs are visually obvious have Wc < 3 A, while objects with Wc > 4 A 
show no sign of these features. W^ib (where "wlb" stands for "weak line blend" ) is centered 
on H9, but has a large value (similar to the Call K line, Wk) even in the absence of H9 
because of the blend of many weak metal lines. Detailed discussions on the properties of 
these indices can be found in Bica etal. (1994) and Storchi-Bergmann etal. (2000). 

Table 4. Window definitions for the Pseudo Continuum 



Pivot A [A] 


Interval [A] 


4 


bx 


3660 


3657-3667 


0.0165 


1.0069 


3780 


3775-3785 


0.0130 


1.0274 


4020 


4009-4020 


0.0161 


1.0019 


4510 


4501-4512 


0.0820 


0.9337 


4630 


4609-4627 


0.1072 


0.9037 


5313 


5307-5317 


-0.0437 


1.0507 



^The coefficient ax is in units of the median flux 
in the 4789-4839 A interval used to normalize the 
spectra, bx is adimensional 
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In Table 5 we present the results for our nuclear spectra. The continuum colors 
PC3Q60/ and -PC4510/-PC402O) corrected for Galactic extinction, are also listed. The 
uncertainties in all these spectral indices were estimated by means of Monte Carlo simula- 
tions. Each spectrum was perturbed 1000 times with Gaussian noise with amplitude equal 
to the rms fluctuation in the 4010-4060 A interval. All indices were measured for each real- 
ization of the noise, and the la uncertainties were computed from the dispersion among the 
1000 perturbed spectra. The median uncertainties are between 0.3 and 0.5 A for all Wx^s 
and 0.03 for the PCsqgq/ PC4020 and -PC4510/PC4020 colors. The largest uncertainties are not 
more than twice these values. 

In Figure 11 we illustrate the relation between the rj = Y, I, I/O and O empirical 
stellar population classes defined in §4.1 and the equivalents widths Wc, Wk, Wq, Wmq 
and the -PC4510/-PC4020 color. The correspondence is excellent, as can be appreciated by the 
locations of different symbols in all panels in this figure. Nuclei with young or intermediate 
age populations are confined to the low W\ and blue color bottom-left corners of Figures 
lla-d, while nuclei dominated by old stars populate the high Wa, red region, with 7] = I/O 
objects in between. Approximate dividing lines can be placed at Wc — 3.5, Wk — 15 
Wg = 9,WMg = 9A and PC4510/PC4020 = 1-55. 

Any of the Wx^s in Figure 11 can therefore be used as a continuous stellar population 
tracer which replaces our discrete /, I/O. O classification. In our previous study of the 
stellar populations in type 2 Scyfcrts (Cid Fcrnandcs ctal. 2001a) Wk was used in this way. 
Of all absorption features measured here, Wk is the strongest one, which makes it also the 
more precise in relative terms. Although the K line is actually produced by late type stars, it 
indirectly traces younger populations by the dilution they cause on Wk- For this reason, in 
practice Wk is nearly as good a tracer of HOBLs as a more direct measure of these features, 
hke Wc (Figure 11a). 

5.2. Spectral Properties in HFS's system 

HFS97 used a set of spectral indices (defined by Whorthey 1992 and Held & Mould 
1994) which also allow a characterization of stellar populations. Seven of the nine absorption 

Fig. 10. — Examples of the pseudo continuum in Bica's system, traced by the automatic 
method developed here (see text). Solid dots mark the pivot continuum points. Vertical 
dotted lines indicate windows of integration of some of the absorption features in this system. 
From left to right: Wc, Wk, Wcn, Wc and Wmq- 
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Table 5. Nuclear Spectral Properties 



Galaxy 


Wc 




Wk 


Wh 


WcN 


Wg 




■^3660 


^4510 


D„(4000) 


H5a 


TT 

W{0 11) 


NGC 0266 


4.6 


18.0 


18.8 


13.5 


15.7 


10.7 


11.6 


0.55 


1.53 


2.03 


-2.1 


5.4 


NGC 0315 


4.3 


17.4 


17.0 


11.7 


15.1 


10.7 


11.3 


0.63 


1.57 


2.00 


-2.8 


8.1 


NGC 0404 


1.6 


9.6 


9.8 


11.3 


7.4 


6.5 


5.3 


0.44 


1.35 


1.50 


4.5 


9.9 


NGC 0410 


4.9 


19.4 


17.6 


12.9 


17.6 


10.6 


12.9 


0.60 


1.55 


2.13 


-3.1 


-1.0 


NGC 0428 


2.0 


10.9 


14.8 


12.6 


8.5 


7.3 


4.4 


0.49 


1.23 


1.65 


2.4 


-1.7 


NGC (J521 


4.5 


18.1 


17.9 


13.4 


17.4 


11.0 


12.3 


0.52 


1.59 


2.01 


-2.7 


-1.3 


NGC 0660 


1.3 


10.8 


14.5 


11.5 


7.3 


7.5 


7.1 


0.47 


1.71 


1.72 


-1.4 


22.4 


NGC 0718 


2.6 


13.1 


13.1 


12.9 


9.5 


8.3 


7.7 


0.52 


1.30 


1.65 


3.1 


1.0 


NGC 0772 


2.5 


10.5 


11.6 


10.0 


11.0 


8.4 


8.2 


0.68 


1.29 


1.46 


-0.3 


0.8 


NGC 0841 


2.7 


13.5 


14.9 


13.0 


10.4 


8.9 


8.5 


0.55 


1.32 


1.73 


1.9 


1.7 


NGC 1052 


4.5 


15.8 


17.3 


8.2 


18.1 


9.9 


11.1 


0.63 


1.61 


1.88 


-1.3 


87.5 


NGC 1161 


5.0 


19.4 


19.0 


13.8 


18.5 


11.6 


12.2 


0.56 


1.63 


2.16 


-2.8 


-1.7 


NGC 1169 


4.3 


16.7 


18.8 


12.8 


15.5 


12.0 


11.0 


0.61 


1.57 


2.08 


-2.4 


7.3 


NGC 2681 


1.8 


11.5 


12.3 


12.8 


8.3 


7.2 


6.8 


0.47 


1.23 


1.57 


4.9 


1.6 


NGC 2685 


5.4 


17.6 


18.7 


13.6 


15.0 


11.6 


10.1 


0.55 


1.53 


2.01 


-1.9 


0.6 


NGC 2911 


4.2 


18.4 


17.7 


13.5 


17.7 


10.9 


11.2 


0.51 


1.68 


2.05 


-2.0 


29.1 


NGC 3166 


3.2 


15.2 


15.9 


13.0 


12.1 


9.4 


8.5 


0.53 


1.39 


1.80 


0.8 


4.7 


NGC 3169 


4.0 


16.6 


17.3 


13.3 


12.9 


9.9 


9.0 


0.50 


1.60 


1.99 


-0.4 


6.2 


NGC 3226 


4.4 


18.2 


18.5 


13.4 


17.6 


11.3 


11.7 


0.58 


1.74 


2.13 


-2.4 


7.7 


NGC 3245 


3.9 


1 5 . 3 


15.2 


11.8 


12.7 


9.6 


10.6 


0.67 


1.45 


1.79 


-1.5 


1.7 


NCX^ 3627 


1.7 


ll.f) 


11.6 


12.7 


7.t) 


7.0 


6.6 


0.48 


1.29 


1.55 


4.4 


7.8 


NGC 3705 


3.0 


12.9 


14.8 


13.1 


10.1 


8.7 


7.6 


0.49 


1.40 


1.73 


2.2 


4.0 


NGC 4150 


2.1 


11.2 


12.6 


12.3 


7.7 


7.4 


4.9 


0.51 


1.19 


1.55 


3.6 


2.3 


NGC 4192 


3.0 


15.3 


16.0 


10.9 


11.8 


9.0 


9.5 


0.53 


1.76 


1.93 


0.1 


13.4 


NGC 4438 


4.1 


17.4 


17.9 


12.6 


15.7 


11.4 


9.9 


0.51 


1.80 


2.00 


-1.9 


57.8 


NGC 4569 


0.6 


7.4 


5.0 


9 . 5 


1.5 


2.7 


3.7 


0.60 


0.94 


1.21 


6.4 


1.2 


NGC 4736 


2.6 


13.6 


12.9 


12.3 


9.9 


8.4 


7.6 


0.53 


1.29 


1.67 


2.8 


-0.1 


NGC 4826 


2.7 


15.2 


14.4 


12.7 


13.2 


9.6 


9.1 


0.52 


1.41 


1.76 


1.0 


4.6 


NGC 5005 


2.7 


13.8 


14.6 


12.1 


10.6 


8.3 


6.7 


0.48 


1.46 


1.74 


2.3 


38.1 


NGC 5055 


3.8 


15.1 


14.1 


11.9 


15.1 


8.7 


10.5 


0.54 


1.47 


1.75 


-0.9 


-1.3 


r<\j\_, 5377 


1.8 


10.7 


8.7 


11.0 


8.0 


6.7 


7.8 


0.53 


1.18 


1.45 


3.7 


4.0 


NGC 5678 


1.1 


7.8 


8.7 


10.3 


3.8 


6.0 


5.6 


0.49 


1 .40 


1 .48 


3.6 


-0.1 


NGC 5879 


3.2 


11.7 


16.3 


12.7 


8.8 


8.7 


7.5 


0.48 


1 .33 


1 . 72 


0.8 


3.2 


NGC 5921 


2.3 


11.2 


11.2 


11.6 


6.5 


6.5 


6.2 


0.53 


1.18 


1.51 


4.6 


5.0 


NGC 5970 


4.1 


14.5 


18.4 


13.3 


10.1 


10.2 


8.0 


0.61 


1.40 


1.86 


-0.8 


-2.5 


NGC 5982 


4.6 


18.4 


18.1 


13.7 


16.1 


10.8 


11.0 


0.57 


1.47 


2.05 


-1.9 


-2.5 


NGC 5985 


4.0 


15.5 


18.9 


12.1 


13.2 


10.6 


9.6 


0.59 


1.51 


1.97 


-2.6 


25.8 


NGC 6340 


4.9 


18.1 


20.0 


14.3 


17.1 


12.1 


11.2 


0.49 


1.92 


2.26 


-2.5 


2.4 


NGC 6384 


4.9 


15.2 


18.6 


12.9 


10.3 


10.7 


9.9 


0.66 


1.45 


1 .92 


-2.3 


-2.7 


NGC 6482 


4.7 


18.7 


18.8 


13.4 


17.5 


11.3 


12.8 


0.63 


1.54 


2.11 


-2.7 


4.4 


NGC 6500 


3.8 


14.5 


15.8 


8.7 


14.5 


10.0 


10.6 


0.74 


1.42 


1.74 


-4.3 


98.2 


NGC 6501 


4.2 


17.7 


16.8 


12.9 


17.7 


10.8 


11.9 


0.56 


1.56 


2.03 


-2.2 


-3.4 


NGC 6503 


1.3 


7.9 


9.5 


10.4 


4.3 


5.8 


6.2 


0.67 


1.15 


1.36 


2.9 


4.5 


NGC 6702 


4.5 


18.5 


18.1 


13.3 


14.0 


10.8 


11.0 


0.61 


1.49 


2.04 


-1.7 


-3.0 


NGC 6703 


4.8 


18.7 


18.5 


13.8 


16.3 


11.5 


11.5 


0.56 


1.60 


2.11 


-2.4 


-1.9 


NGC 6951 


3.7 


14.1 


16.4 


11.7 


15.7 


11.0 


10.2 


0.52 


1.43 


1.70 


-0.2 


21.3 


NGC 7177 


3..') 


14.5 


16.6 


12.5 


11.7 


10.3 


9.2 


0.58 


1.52 


1.80 


-1.3 


11.2 


NGC 7217 


5.0 


18.7 


19.2 


14.0 


17.2 


11.5 


12.0 


0.53 


1.68 


2.07 


-2.6 


12.2 


NGC 7331 


4.2 


16.8 


18.0 


13.6 


15.3 


11.0 


10.9 


0.54 


1.55 


1.99 


-1.0 


-1.0 


NGC 7626 


4.9 


19.2 


18.1 


13.9 


17.9 


11.0 


12.6 


0.58 


1.58 


2.12 


-3.3 


-1.2 


NGC 7742 


4.1 


15.9 


17.1 


13.1 


13.1 


10.6 


9.6 


0.52 


1.50 


1.90 


-0.9 


2.2 


NGC 3367 


0.3 


1.8 


2.6 


1.1 


1.5 


1.4 


1.6 


1.07 


0.86 


1.05 


-3.0 


13.3 


NGC 6217 


0.3 


3.9 


2.5 


5.3 


2.7 


2.2 


2.3 


0.75 


0.88 


1.13 


3.1 


5.9 


NGC 0205 


1.2 


8.4 


7.1 


11.5 


4.3 


5.2 


3.3 


0.51 


1.07 


1.37 


6.3 


-2.0 


NGC 0221 


4.0 


15.5 


17.9 


13.9 


11.9 


10.6 


8.1 


0.60 


1.41 


1.86 


-0.9 


-3.4 


NGC 0224 


4.8 


18.6 


17.5 


13.7 


18.5 


10.3 


12.0 


0.63 


1.54 


2.06 


-2.9 


-4.1 


NGC 0628 


2.9 


13.1 


16.1 


12.7 


8.5 


8.5 


6.8 


0.64 


1.35 


1.74 


0.5 


-2.4 


NGC 1023 


5.4 


20.8 


19.4 


14.5 


19.1 


11.9 


12.3 


0.57 


1.63 


2.26 


-3.7 


-4.6 


NGC 2950 


4.6 


18.7 


17.4 


13.2 


17.0 


10.6 


11.3 


0.59 


1.48 


2.04 


-2.0 


-3.0 


NGC 6654 


4.7 


16.4 


18.5 


13.4 


16.6 


12.5 


12.1 


0.61 


1.58 


2.09 


-1.6 


-5.4 



Note. — Col. (1): Galaxy name; Cols. (2—10): Equivalent widths of seven absorption features and two colors, all in Bica's system. Col. 
(11): 4000 A break index (Balogh et al. 1999). Col. (12): H^a equivalent width of Worthey & Ottaviani (1997). Col. (13): [Oil] equivalent 
width of Balogh etal. (1999). All equivalent widths are given in A 
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equivalent widths and one of the two colors measured by HFS97 are in the wavelength range 
covered by our data: the G-band, H7, Fe4383, Ca4455, Fe4531, Fe4668, Ef3 and c(44 - 49). 
We have measured these indices in our spectra following the recipes in HFS97. This was 
done mainly for comparison purposes. 

In Figure 12a we compare the equivalent width of the G-band for our HFS-aperture 
spectra with the measurements of HFS97. The mean difference is just 0.03 A, and the rms 
is 0.6 A, within the measurement errors of both studies. Similarly, Figure 12b compares the 
values of c(44 — 49) in these two studies. The difference in c(44 — 49) is just 0.03 ± 0.04, 
again within the errors. 

In Figure 12c we compare the equivalent widths of the G band as measured in the HFS 
system and Bica's Wg, while in Figure 12d we plot the c(44 — 49) color against Bica's Wk- 
All quantities are measured from our nuclear spectra. The relatively tight relations between 
these independent measurements reinforces the impression from Figure 11 that different 
stellar indices in galaxy spectra are highly correlated. A corollary of these correlations is 
that a single index (say, Wk) is enough to provide a good first order characterization of 
the stellar populations in our sample galaxies. This relatively high degree of redundancy in 
galaxy spectra is also detected by studies of principal component analysis (e.g., Schmidt et al. 
1991; Sodre & Cucvas 1997), which show that a large fraction of the variance in galactic 
spectra is described by a single principal component. The physical property behind such 
an approximately mono-parametric behavior is the age of the dominant stellar population 
(Schmidt etal. 1991; Ronen, Aragon-Salamanca & Lahav 1999). 

The strong inter-relations between different stellar population indices can be used to 
extend results drawn from our sample to the full HFS sample. For instance, any correlation 
found between Wc or Wk and emission line properties for our galaxies will also appear 
with the 14^(G-band) HFS index, which has been measured by HFS97 for most of their 162 
LLAGN. Analogously, statistics related to our Y, I, I/O and O stellar population classes 
can be readily extended to the full HFS sample exploring the relation between rj and, say, 
ly(G-band). 

Fig. 11. — Relations between the equivalent widths of four absorption lines and one color 
in Bica's system. Empirical stellar population classes Y, /, I/O and O are represented 
by different symbols. Squares, triangles, circles and stars represent rj — O, I /O, I and Y 
respectively. (The two stars correspond to NGC 3367 and NGC 6217, two Starburst nuclei 
in our comparison sample.) Dotted lines approximately separate classes rj — Y and / from 
rj — I/O and O. Mean error bars are indicated in the bottom right corner of the figures. 
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5.3. Other spectral indices: D„(4000), USa and W(0 II) 

Besides the Call K line and the HOBLs, the interval between 3500 and 4200 A contains 
other features of interest. The 4000 A break and ¥16, for instance, have been recently used as 
stellar population diagnostics for SDSS spectra of both normal and active galaxies (Kauffman 
etal. 2002, 2003). We have measured these two indices for our spectra following the Balogh 
etal. (1999) definition of Dn(4000) and the Worthey & Ottaviani (1997) definition of RSa- 
These are the same indices used in the SDSS work. We have not corrected RSa for emission 
as this is insignificant for most of our LLAGN. We further measured the [Oil] equivalent 
width 14^(0 II) according to the recipe of Balogh etal. (1999). All spectra were corrected for 
Galactic reddening prior to these measurements. 

These indices are listed in the last three columns of Table 5. Note that RSa > 
corresponds to absorption while W{0 II) > corresponds to emission. Note also that 

negative W{0 II) is obtained in many cases, but this does not necessarily means a non- 
detection of [Oil]. For instance, [Oil] is clearly seen in the starlight subtracted spectrum 
of NGC 7626, while the W{0 II) index, which is measured over the total spectrum, yields 
—1.2 A. As for other emission lines in LLAGN, the accurate measurement of [Oil] requires 
a careful starlight subtraction. 

In Figure 13a we examine the relation between D„(4000) and USa- The global distri- 
bution of LLAGN in this diagram is similar to that obtained for SDSS galaxies (Kaufmman 
etal. 2002). Interestingly, our rj — I LLAGN, which are nearly all weak-[OI] sources, lie 
above the locus of steady star formation models of Kaufmann et al. (2002), in a region occu- 
pied by galaxies which have experienced bursts of star formation in the past ~ 10^ yr. This 
is consistent with the strong HOBLs found in these systems and the considerations of §4.5. 
The location oi r] = I weak- [01] sources in Figure 13a is also similar to that of AGN with 
L[oiii] < 10^ Lq and TO-like emission line ratios in Kauffman etal. (2003). Notwithstanding 
these similarities, it is important to observe that while the SDSS work pertains to scales of 
several kpc, our data traces stellar populations on scales of typically 70 pc. 

Figure 13b illustrates the strong relation between Bica's Wk and £)n(4000), which is 

Fig. 12. — Top: (a) Comparison of the equivalent width of the G band and (b) the c(44 — 49) 
color for our HFS-aperture spectra and the measurements of HFS97. Both indices were 
computed according to the definitions in HFS97. Diagonal fines indicate y — x and the ±2 

sigma range. Bottom: Relation between Wq and Wk, measured through Bica's system, and 
two HFS indices: (c) the G-band equivalent width and (d) c(44 — 49). Symbols as in Figure 
11. 
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not surprising given that the Wk window is contained in blue window of D„(4000). The 
dashed hne indicates the quadratic fit: 

L>„(4000) = (1.013 ± 0.006) + (0.031 ± 0.001) W^x + (0.00124 ± 0.00007) (2) 



6. The Link Between Stelleir Population and Emission Line Properties in 

LLAGN 

Several times in §4 we have noted an apparent connection between stellar population 
features revealed by our spectroscopy and the weak/strong [01] classification. In this section 
we use the spectral properties measured in §5 in conjunction with the [01]/Hq; ratio from 
HFS97 to investigate this connection in a quantitative manner. 

One way to investigate this relation more closely is to correlate [01] /Ha with the Wc 
index, a useful tracer of the presence of HOBLs. The strong tendency of HOBLs to appear 
preferentially in weak- [01] nuclei noted in §4.3 should reveal itself as a correlation between 
[OI]/Hq; and Wc- Similar correlations should exist between [01] /Ho; and the equivalent 
widths of metal absorption lines, since HOBLs come with an associated blue continuum 
which dilutes spectral features like the K hne and the G band. 

Figure 14 confirms these expectations. The figure plots [01] /Ha against our nuclear 
measurements of Wc, Wk, Wc and Wmq (Table 5). What is seen in these plots is not so 
much a correlation in the sense of points scattered around a line, but a marked dichotomy: 
Nearly all nuclei with weak absorption lines are weak- [01] emitters, whereas nuclei with 
strong metal lines span the full range of [01] /Ha. In graphical terms, the top-left regions of 
all panels in Figure 14 are remarkably empty! 

This is not the first time in this paper that [01] /Ha is plotted against a stellar population 
index. In Figure 2 we plotted [01]/Ha versus W{G band) for the full HFS sample. There 
we see the same "inverted L" shape of any of the panels of Figure 14. Overall, these figures 
spell out a clear message: The Narrow Line Region of LLAGN evidently knows what value 
of [01] /Ha it can or cannot assume for a given stellar W\. This finding implies a strong 
connection between stellar populations and gas excitation processes in LLAGN. 

Fig. 13. — Relations between D„(4000), B.Sa and Wk- Symbols as in Figure 11. The dashed 
line in (b) shows a second order polynomial fit. 
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As suggested by the referee, an alternative reading of this result could be that galaxies 
with important circumnuclear star formation have an excess of Ha emission in the central 
few arcsec (note that the [01] /Ha values are from the 2" x 4" measurements by HFS97), 
resulting in a lower [01] /Ha ratio. We have studied the aperture dependence of the Ha 
flux in order to check whether the link between [01] /Ha and the stellar population is not a 
consequence of an excess of Ha flux in the larger HFS97 apertures, compared to the smaller 
apertures (1" x 1.1") used in this work to derive the stellar population properties. We have 
retrieved from the HST archive STIS (G750M or G750L) spectra of 27 out of the 28 LLAGN 
analysed in Paper II. We find that Ha is very compact in all strong- [01] LLAGN. In weak- 
[01] LLAGN, the geometry of Ha is more diverse: in 7/22 Ha is compact, and in 14/22 more 
than 50% of the flux (measured through a 0.2" x 4" extraction) is inside the central 0.2" x 1". 
On the other hand, we find that the extension of the Ha emission does not correlate with the 
type of nuclear stellar population. In addition, we find that there is no correlation between 
the Ha luminosity and the [01] /Ha ratio (both measured by HFS97). We thus conclude that 
the lower values of [01] /Ha in weak- [01] LLAGNs are not due to an excess of Ha emission 
provided by recent star formation in spatial scales of a few arcsecs. 

A compelling visualization of the relation between [01] /Ha and the stellar population 
properties is presented in Figure 15. In this figure we divide objects into bins in stellar Wx 
and within each bin we compute the fraction of objects with [OI]/Ha < 0.17 (filled areas), 
0.17 < [01] /Ha < 0.25 (hatched areas) and [01] /Ha > 0.25 (empty areas). The panels with 
Wk and Wc are for our sample, while the panel with l^(G-band) uses data tabulated in 
HFS97 for all their LLAGN. 

The predominance of weak [01] emitters among objects with low stellar absorption 
equivalent widths is obvious to the eye. Of the 19 LLAGN with Wk < 15 A, 17 have 
[OI]/Ha < 0.25. Similarly 18 of the 21 nuclei with Wc < 3.5 A have [OI]/Ha < 0.25. 
Hence, ~ 90% of objects with weak absorption fines are also weak [01] emitters. Note that 
the above limits in Wk and Wc correspond approximately to the border line between our 
f] = I and 1] = I/O or O stellar population classes (Figure 11), so nearly identical statistics 
are obtained comparing 1] with [01] /Ha. Likewise, since rj = I systems are also those with 
HOBLs, one can alternatively express the above statistics in terms of a [01]/Ha-H0BLs 

Fig. 14. — Comparison of absorption equivalent widths, which trace stellar populations, with 
the [01] /Ha ratio, an AGN indicator. The [01] /Ha = 0.25 border-line which separates weak 
from strong- [01] nuclei is indicated. The vertical dotted lines arc the same as in Figure 
11; they separate objects with a strong intermediate age population {rj = I) from objects 
dominated by old stars {r] — O). Symbols as in Figure 11. 
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connection. 

For the full HFS sample, - 80% of the W{G band) < 4.5 A sources have [OI]/Hq; < 
0.25 (Figure 15c). This same fraction results if one uses the HFS97 measurements for the 
51 LLAGN in our sample, which reinforces our conclusion that there are no significant 
differential selection effects between the two samples. 

It is important to draw attention to the fact that although most objects with HOBLs 
and diluted metal lines are weak in [01], the converse is not true. There is a numerous 
population of weak- [01] sources without conspicuous HOBLs and with stellar population 
indices typical of old, ~ 10 Gyr populations. Among objects with Wk > 15 A, for instance, 
~ 50% have [OI]/Ha < 0.25. These are the objects located at the bottom-right corners of 
the "inverted L" distributions in Figures 2 and 14, also represented by the filled and hatched 
areas in the large Wx bins in Figure 15. 

6.1. Interpretation 

The shape of the distribution of points in Figure 14b lends itself to division into four 
regions, delimited by [OI]/Hq; ~ 0.25 in the vertical (emission line) axis and Wk ~ 15 A 
in the horizontal (stellar population) axis. Similar divisions can be placed at Wc ~ 3.5, 
Wg ~ 9 and WMg ~ 9 A in the other panels of Figure 14, and at 14^(G-band) ~ 4.5 A in 
Figure 2. The top- left of these four quadrants is essentially empty. What kind of objects 
populate the remaining three quadrants? 

Objects with clear signs of stellar populations of 10^ yr or less (ie, those with HOBLs, 
diluted metal lines and blue colors) all live in the bottom-left quadrant. They arc all wcak- 
[01] 's according to our definition. A plausible reading of this fact is that in these objects 
stars are directly linked to the gas excitation by some as yet undefined mechanism (see §1 
for a menu of possibilities), resulting in higher Ha per [01] photon ratio than in strong- [01] 
LINERs. These sources could be called "Young-TOs'\ where "young" actually refers to 
populations of 10^ yr of less. Note that the stellar Wa's in these sources are intermediate 
between those of starburst and early type galaxies. In other words, they are transition 

Fig. 15. — Fraction of objects in each of three ranges in [01] /Ha as a function of different 
stellar absorption equivalent widths. Filled areas are used for [01] /Ha < 0.17, hatched areas 
for 0.17 < [OI]/Ha < 0.25 and empty boxes for [OI]/Ha; > 0.25. Numbers at the top indicate 
the number of objects in each PFx-bin. The left and middle panels are for the present sample 
of 51 LLAGN, while the right panel is for the full HFS sample. 
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objects also in terms of their stellar populations. 

In the opposite end of the [OI]/Hq; versus Wx diagrams, stars cannot play an important 
role in the gas excitation, since photoionization by hot stars, be they young (ie, a starburst) 
or old (post-AGB stars), produces smaller values of [OI]/Hq; than observed in AGN. The 
absence of optically significant < 10^ yr populations corroborates this interpretation. Pho- 
toionization by an AGN is likely the dominant emission line mechanism in these sources. 
Variations in nebular conditions or in the ionizing spectrum can be invoked to explain the 
range of [01] /Ha values spanned by the data. A fitting denomination for nuclei in this 
top-right quadrant is "Old-LINERs". 

The bottom right quadrant of Figure 14 is populated by objects with TO-like emis- 
sion line spectra but no strong optical signs of young or intermediate age stellar popula- 
tions. These "Old-TOs" may well belong to the Old-LINER family, in which case their weak 
[01] /Ha should perhaps be explained by an ionizing radiation field softer than for strong- 
[01] LINERs. Alternatively, they may be simply weak Young- TOs whose < 10^ yr stars are 
completely overwhelmed by the older population. This is hkely the case of NGC 5055. The 
young starburst which clearly dominates the UV hght of this galactic nucleus (Maoz et al. 
1998) leaves no trace in the optical (Table 3), although our starhght decomposition indicates 
that 21% of the 4020 A flux originates in an intermediate age population. This non-ncgligiblc 
contribution, however, causes little dilution of the metal lines. The lower right quadrant of 
Figure 14 probably contains other such objects. 

To summarize, our results suggest a scenario where stars younger than ~ 1 Gyr somehow 
participate in the gas ionization of Young-TOs, whereas in Old-LINERs an AGN is the 
dominant ionizing agent. The situation is less clear for Old-TOs, which may comprise a 
mixture of the above objects. Interestingly, "Young-LINERs" do not seem to exist. AGN 
taxonomy surely does not need more sub-classes. Yet, LLAGN have long been known to 
comprise a heterogeneous class, and, as discussed in §1, models involving AGN- and stellar- 
driven ionization can both explain their emission line properties. This empirically motivated 
scenario comes in the sense of trying to elucidate which of these two broad physical processes, 
AGN or stars, best fits a given source. It is equally important to clarify that we do not claim 
Young-TOs to be AGN-lcss. After-all, the current common wisdom is that all galaxies with 
bulges harbor an AGN, be it dead, alive or somewhere in between. In Young-TOs the AGN 
contribution to the gas ionization would be small. Conversely, ionization by stellar processes 
in Old LINERs would be overwhelmed by the AGN. 

A potential caveat in identifying nuclei with weak metallic lines with stellar processes 
is that dilution by a quasar-like non-stellar continuum can also produce the same effect. 
Presumably, this should be more important for type 1 LLAGN, ie, LI. 9 and T1.9 in the 
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notation of HFS97. This effect indeed happens in HST spectra of some objects (Paper II). 
Interestingly, the two objects with W{G band) between 3.5 and 4 A and [OI]/Ha ~ 0.6 in 
Figure 2 which intrude shghtly into the "zone of avoidance" arc NGC 841 and NGC 5005, 
both classified as LI. 9. Although it is conceivable that this effect operates in a few cases, 
wc find that type 1 and type 2 sources are well mixed in all our diagrams, so this cannot 
be a major effect in our ground based observations. Furthermore, we point out that HOBLs 
were detected in both NGC 841 and NGC 5005, as well as in NGC 2681, another L1.9 with 
relatively weak metal absorption lines. It therefore seems more likely that the main diluting 
agent even in these type 1 LLAGN is the continuum associated with stellar populations of 
< 10^ yr rather than with a quasar- like continuum. 



6.2. Evolution and Contrast 

As time passes, HOBLs and the continuum of < 10*^ yr populations gradually disappears, 
and metal absorption features become progressively stronger. Passive evolution thus moves 
nuclei in left to right direction in Figures 2 and 14 on time scales of a few Gyr. In the notation 
of the previous section, objects in the left half of Figures 2 and 14 are nearly all Young TOs. 
If [01] /Ha remains constant during this period, stellar evolution would gradually turn a 
Young-TO onto an Old-TO. This would correspond to a horizontal evolutionary track onto 
any of our [01] /Ho; versus Wx diagrams. 

However, the hypothesis of constant [OI]/Hq; does not seem warranted by the "inverted 
L" distribution of LLAGN in Figures 2 and 14, which suggests that the presence of < 10^ 
yr populations affects the [01] /Ha ratio. For didactic purposes, lets consider a situation 
in which an AGN is entirely responsible for the production of [01] photons, while Ha is 
produced both by the AGN and < 10^ yr stars. In this case the stellar contribution to 
[01] /Ha acts just in the denominator of this ratio. Evolution would thus move objects both 
to the right and upwards along curved mixing lines in Figures 2 and 14 as the stellar powered 
Ha emission fades. Considering that the AGN part of [01] /Ha can assume a range of values 
(due, eg, to different nebular conditions), one arrives at a scenario where stellar evolution 
takes points in the bottom-left to any possible height in the right part of these figures. In 
other words, Young-TOs can end up as either Old-TOs or Old-LINERs. 

It thus seem plausible that evolution plays a role in defining the looks of LLAGN spectra, 
both in their stellar and emission lines properties. Physical models for the evolution of 
[01] /Ha must bear in mind that the intensity of [01] is not only driven by the ionizing 
radiation field; it also depends strongly of the nebular geometry. Thus, variations of the 
ionization parameter and/or the density structure would also affect the [OI]/Ha ratio. Such 
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variations may be tied to the evolution of stellar populations. In starbursts, for instance, 
[OI]/Hq; is expected to increase with age as the bubble produced by the action of stellar 
winds and SNe expands (Stasinska, Schaerer & Leitherer 2001; Stasinska & Izotov 2003). 

Besides evolution, contrast effects must also be present, as can be illustrated by the 
case of NGC 6500. Judging by the relative contributions of the Y, I and I/O components 
in our starlight decomposition (25, and 4% of the 4020 A flux, respectively), this is the 
youngest LLAGN in our sample. Yet, since 71% of the optical light comes from old stars, 
this galaxy does not quite make it into our observational definition of Young-TO, whereas 
from the above discussion of evolutionary effects one could naively expect it to be located at 
the bottom left of our [01] /Ha versus W\ diagrams. Further examples of contrast effects are 
NGC 5055 and NGG 404, whose optical spectra reveal no sign of the young massive stars 
unambiguously detected at UV wavelengths by Maoz et al. (1998). This starburst component 
is simply buried under the much stronger contribution of intermediate age and older stars. 

A quantitative assessment of these ideas is beyond the scope of this paper. Future 
communications will explore these issues by means of stellar population synthesis, spectral 
energy distribution modeling and comparisons to Seyfert and Starburst galaxies. 

7. Conclusions 

In this paper we presented the first results of a survey of the stellar populations of low 
luminosity AGN. Nuclear spectra covering the inner ~ 70 pc in radius were presented for 51 
LINERs and LINER/HII Transition Objects plus a comparison sample 9 non-active galaxies. 
The stellar population properties of LLAGN were quantified by means of a comparison with 
non-active galaxies as well as by measurement of a suite of spectral indices over the 3500-5500 
A interval. 

Our main results can be summarized as follows: 

1. Young, < 10^ yr starbursts contribute very little to the optical spectrum of LLAGN. 
Even in the few nuclei known from UV observations to contain young massive stars, 
these populations participate with no more than 30% of the optical continuum light, 
and sometimes are altogether undetected due to a much stronger contribution from 
other stellar populations. 

2. Intermediate age, 10^-10^ yr populations, on the other hand, are found in at least 1/3 
of all LLAGN. These populations are easily characterized by pronounced HI Balmer 
lines in absorption, diluted metal lines and relatively blue colors. 
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3. We found a very strong empirical connection between the stellar populations and emis- 
sion line properties of LLAGN. Nearly all nuclei containing substantial populations of 
young and/or intermediate age stars are weak [01] emitters, with [OI]/Hq; < 0.25, 
whereas only ~ 10% of the sources with stronger [01] contain such populations. The 
mean stellar age of weak- [01] LLAGN is thus several Gyr younger than that of strong- 
[01] LLAGN. Nuclei dominated by older stars, however, span the full range of [OI]/Ha 
values, from TOs to LINERs. 

4. Given a few Gyr, TOs with conspicuous intermediate age populations will naturally 
evolve to either "Old-TOs" or "Old-LINERs" . 

Perhaps the major lesson from this study is that the information contained in the stellar 
populations of active galaxies can provide valuable clues on the nature of these systems. 
Combining this information with emission line data, we were able to establish an inequivocal 
connection between stellar populations and the gas excitation mechanism in LLAGN. The 
physical origin of this connection, however, has not been unveiled in this work. By analogy 
with previous investigations (particularly of type 2 Seyferts) , starburst activity is the prime 
contender to explain the relation between [01] /Ho; and stellar population properties reported 
in this work. Direct evidence for young stars in LLAGN is nevertheless too scant to fully 
subscribe this interpretation at this stage, specially because other scenarios (such as SN 
shocks and photoionization by evolved post-starburst populations) may potentially explain 
our findings. Future papers in this series will address these and other related issues exploring 
complementary spectroscopic and imaging data, as well as detailed modeling of the stellar 
populations in these sources. 
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